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Abstract A stochastic background of primordial gravitational waves may substantially 
contribute, via the Sachs-Wolfe effect, to the large-scale Cosmic Microwave Background 
(CMB) anisotropies recently detected by COBE. This implies a bias in any resulting deter- 
mination of the primordial amplitude of density fluctuations. We consider the constraints 
imposed on n < 1 ("tilted") power-law fluctuation spectra, taking into account the con- 
tribution from both scalar and tensor waves, as predicted by power-law inflation. The 
gravitational-wave contribution to CMB anisotropies generally reduces the required rms 
level of mass fluctuation, thereby increasing the linear bias parameter, even in models where 
the spectral index is close to the Harrison-Zel'dovich value n = 1. This "gravitational- 
wave bias" helps to reconcile the predictions of CDM models with observations on pairwise 
galaxy velocity dispersion on small scales. 

Subject headings: cosmic background radiation — cosmology — early universe — galaxies: 
formation 
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The recent detection of large angular scale CMB anisotropies by the COBE satellite 
(Smoot et al. 1992) opens a window to the understanding of the physics of the early uni- 
verse: in particular, it provides strong constraints on models for the origin of primordial 
perturbations. Inflation is probably the simplest and most motivated of such models: per- 
turbations are generated in a causal way by zero-point quantum fluctuations which are then 
magnified by the accelerated universe expansion to cosmologically observable scales. The 
determination of the rms fluctuation amplitude consistent with the COBE measurements 
in the frame of various galaxy formation scenarios (e.g. Wright et al. 1992; Efstathiou, 
Bond & White, 1992; Schaefer & Shafi 1992) has however shown that a quite high fluctu- 
ation level is required, which, in the standard Cold Dark Matter (CDM) scenario causes 
excessive small-scale power. In fact, a relevant quantity for all galaxy formation scenarios 
is the linear bias parameter, defined as the inverse of the rms linear mass fluctuation on a 
sphere of 8 h~ x Mpc (h is the Hubble constant in units of 100 km sec -1 Mpc -1 ; we take 
h = 0.5): the COBE results imply b 0.8, for the scale-invariant (n = 1) CDM case. Low 
bias (i.e. more evolved) CDM models lead to better agreement with observations on large- 
scale flows (e.g. Bertschinger et al. 1990), but imply an excess of velocity dispersion on 
small scales when compared to observations on pairwise galaxy velocity dispersion in the 
CfA redshift survey (Davis & Peebles 1983). Also, the slope of the galaxy two-point func- 
tion, determined in numerical simulations, becomes too steep. Both of these drawbacks 
can be alleviated by resorting to a velocity bias (e.g. Couchman & Carlberg 1992). 

The above determination of b is however only valid if the large angular scale tem- 
perature anisotropies, detected by COBE, are totally due to density perturbations (scalar 
modes), which perturb the last-scattering surface via the Sachs- Wolfe effect (Sachs & 
Wolfe 1967). However, a stochastic background of primordial gravitational waves (tensor 
modes) originated during inflation also contributes to this effect. A rough estimate of 
the anisotropies originated by scalar perturbations is (^-)s ~ ^H5(p/(p : where (p is the 
inflaton field, 6<p its fluctuation and H the Hubble constant during inflation; these quan- 
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tities have to be evaluated at the time when the scales relevant to the large-scale CMB 
fluctuations crossed the Hubble-radius, i.e. about 60 e-foldings before the end of the in- 
flationary expansion. The anisotropy originated by tensor perturbations is (4j£)r ~ f 8(p, 
with k = \^8nG. Thus, we can easily obtain an approximate measure of their relative 
contribution by {™) s / {%f)t ~ \H/k<p\hc- During inflation H 2 = ^-(V((p) + ±v? 2 ), 
where V(cp) is the effective inflaton potential. In many models, such as chaotic (Linde 
1983) or new inflation (Linde 1982; Albrecht & Steinhardt 1982), a slow-rollover approx- 
imation holds, (p 2 <C V((p) ks 3H 2 /k 2 ; in these cases the contribution of tensor modes to 
5T/T is much smaller than that due to scalar ones. However, in other models, such as 
power-law inflation (Abbott & Wise 1984a; Lucchin & Matarrese 1985), a slow-rollover 
approximation is not necessarily required and one can have (p 2 ~ V(cp). More in general, 
the minimal requirement on the inflaton dynamics is that it should lead to accelerated 
universe expansion, which implies ip 2 /V(ip) = e < 1 and (^-)s/(^-)t ~ \\J\ + ^ > "^r- 
It should be pointed out, however, that the possibility to ascribe most of the 8T/T signal 
to gravitational waves (Krauss & White 1992) is restricted to models where the fluctuation 
spectrum is non-scale-invariant. 

To be more specific, let us consider the power-law inflation case, which has the 
advantage of being fully analytically tractable. In this case the universe expansion factor 
reads a(t) = a*[l + (H*/p)(t — t+)] p ~ t p (where a* and H+ refer to an arbitrary time 
t+ during inflation), with p > 1, and the inflaton field is assumed to have an exponential 
potential (Lucchin & Matarrese 1985), V((p) oc exp(— Xnip), with < A = \f2/p < y/2. 
Such an exponential potential also describes the dynamics of extended inflation models (e.g. 
La & Steinhardt 1989; Kolb, Salopek & Turner 1990). Moreover, for any inflation model 
where one scalar field rolls down a smooth potential, the evolution during the small range 
of e-foldings relevant for large-scale CMB anisotropies can be approximated by a power- 
law; thus, our results have a quite general validity. In such a case, the resulting power- 
spectrum of density perturbations at Hubble-radius crossing is proportional to k 2a ~ s , 

3 



with a. = 1/(1 — p) < 0. A stochastic background of gravitational waves is also produced, 
with the same spectral behaviour at Hubble-radius crossing. In what follows we shall 
parametrize both of these spectra by the index n = 2a + 1 = (p — 3)/(p — 1) < 1, which 
for scalar modes (but not for tensor ones!) gives the spectral slope on constant time 
hypersurfaces before recombination, i.e. the so-called primordial spectral index. Note 
that, the limit p — > oo, corresponding to the de Sitter case, gives a — > 0, or n — > 1, i.e. the 
Harrison-Zel'dovich fluctuation spectrum. For a general value of p (or n) one can obtain 
the estimate (^-)s/(^-)t ~ \\J 2 ^i-^ ~ = "^F' thus, for values of n not too far from 
unity, as required by the COBE results, we conclude that gravitational waves may make a 
significant contribution to large angular scale CMB anisotropies. 

Let us now provide a more detailed analysis of the problem. We can perform 
the usual expansion of temperature fluctuations in spherical harmonics ^-(9,<p) = 
m a£ m Yi rn (9, (j)), where the multipole coefficients take independent contributions from 
both scalar and tensor modes: a^ m = as,£ m + a T,tm- Note that, even though we wrote 
both the scalar and tensor modes as being proportional to the same field fluctuation Sep, 
they actually refer to independent quantum field fluctuations, namely the inflaton and one 
polarization state of the graviton, which simply have the same rms value. The squared 
multipole amplitudes a\ = J2 m \ a em\ 2 have expectation values (af) = (df)s + (o>f)T, com- 
ing from both scalar and tensor perturbations. In the simplest case that both the inflaton 
and the graviton fluctuations have random phases, the multipoles ag are Rayleigh dis- 
tributed in 21 + 1 "dimensions" (Abbott & Wise 1984b; Fabbri, Lucchin & Matarrese 
1987), with cosmic variance ^2/(21+ Tj(aj). The result for the scalar case is, in a flat 
universe, (aj) s = ^rp 1 f A*(k)j$(2k/H ) (e.g. Bond & Efstathiou 1987), where 
A$(fc) = T[-^V<f>{k) is the power per logarithmic wavenumber of the peculiar gravitational 
potential $ and V$(k) its power-spectrum. In the power-law inflation case this relation 
can be analytically integrated (e.g. Fabbri, Lucchin & Matarrese 1987; Lyth & Stewart 
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1992) to give 



2 (27r) 4 3-n (H \ n 1 C(n)T(3-n)T(£+(n-l)/2) 



25 1-n v V 2 / 2 2 - n r 2 (2-n/2)r(£- (n-5)/2)' 

where the factor C(n) = 7r42Tt 1 + i a 2 ( ^fz~ ) 3-T T 2 (2 — n/2) is related to the power-spectrum 
of the inflaton by V v {k) = 4nC{n)k n ~ l . The constants a* and H± could be easily related to 
physical observables by matching the inflationary kinematics to the subsequent radiation 
and matter dominated eras. For the gravitational-wave contribution the result is (e.g. 
Abbott & Wise 1984a; Fabbri, Lucchin & Matarrese 1986) 

(a^ = 144 7r 5 G(2£ + l)|±||^ n l C{n) dkk- 2 lj{k), (2) 

where 

.... f" , Je+i/2{k-y)J 5/2 (y) 

and rjE and rjo are the conformal time at the recombination and at the present epoch. In 
this case, the integration must be numerically performed. However, in the range of values 
of n of interest for us, there is a nice property of tensor multipoles (e.g. Fabbri, Lucchin 
& Matarrese 1987) which makes it possible to relate them to the scalar ones in a simple 
manner. By numerically integrating (af)T in Eq.(2) one can show that the ratio Dg(n) = 
{aj)T/(aj)s: for £ > 2, is independent of £ to a very good approximation, Dg(n) D(n). 
From the plots of Figure 1 this can be seen to hold in the spectral range 0.5 ^ n < 1, with 
better than 10% accuracy. Figure 1 also shows that, for n ^ 0.8, gravitational waves give 
the main contribution to the multipoles, while density perturbations dominate for larger 
n values. The ratio of the tensor to the scalar contribution to the quadrupole is larger 
than that due to the higher order multipoles. In the high £ limit (£ 1) it is also possible 
to obtain an approximate asymptotic form for the tensor multipoles, along the lines of 
the computation by Starobinskii (1985) for n = 1. Using Eq.(12) of Starobinskii as an 
approximation for lf(k) we obtain 



{aj) T * 288^C(n)F(n)[^ U+-) . (3) 



with 
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dx 



l-^ 2 -^(l-x 2 )o:ln ^ +^r(l-^ 2 ) 
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which shows the same asymptotic dependence on £ as the scalar multipoles. The numerical 
value of the ratio of the tensor to the scalar components is 5% larger than those plotted 
in Figure 1 for £ = 30, for values of n > 0.5, and 10% larger, for values of n between 
0.5 and 0. We can then write any £ > 2 multipole in terms of the scalar contribution 
only, (a 2 ) = (a 2 )s(l + D(n)). This makes it easier to obtain the bounds imposed by the 
COBE determination of the angular correlation function C{9) of temperature fluctuations 
to the complete scalar + tensor case. After dipole and quadrupole subtraction, one has 
C(6) = (1 + D(n)) £ z>2 (AT 2 ) s W 2 (£)P,(cos#), where ( AT 2 )g = /±7r){aj) s , with T = 
2.735 ± 0.006 K the mean temperature of the CMB radiation (Mather et al. 1990) and 
W(£) = exp[-(l/2)(£(£+l)/17.8 2 )] the appropriate filter function for the DMR experiment 
(Smoot et al. 1992). Thus considering also the gravitational-wave contribution does not 
affect the best-fit on the primordial spectral index n = 1. 15^0^5 (Smoot et al. 1992, 
Wright et al. 1992). This corresponds to a limit on the scale factor expansion power p > 5 



and on the scalar field potential coupling constant A < a/2/5. It affects, instead, the best- 
fit on the amplitude, since the total rms-quadrupole-normalized amplitude Q rms -ps = 
16.3 ± 4.6 //K can now be written as 



Qrms — PS — (Qrms-Ps)sV 1 + D (n)- (5) 

The normalization of primordial density fluctuations required to fit the COBE data 



gets modified by the same factor. In particular one gets b(n) = bo(n)^yi + D(n), where 
bo represents the same quantity calculated disregarding the gravitational-wave contribu- 
tion to CMB fluctuations. Figure 2 shows the ratio b/b as a function of n; notice that 
this "gravitational-wave bias" is independent of the galaxy formation scenario, i.e. of 



the transfer function. A very good fit is b/b ~ \J 1 \}^ n ■ The value of b (n) in the 



frame of the CDM scenario, taking into account the (best-fitted) amplitude and the 
errors given by the COBE results on C(6), can be easily derived by properly scaling 
the fits by Vittorio, Matarrese & Lucchin (1988) [their Eqs.(23) and (24)]: one obtains 
b (n) « 0.75 x io 1285 ( 1 -^)(l ± 0.3). From the COBE detection of the 10° anisotropy, 
a S ky(lO°) = 30 ±5 fiK, we similarly get b (n) « 0.82 x 10 1 - 15 ( 1 " n )(l ± 0.2), in good agree- 
ment with Adams et al. (1992). By taking the weighted (by inverse variance) average of 
the two determinations and considering the gravitational-wave correction we finally get 
the estimate 



V 3 — n 

As it is clear from Figure 1, gravitational waves give an even larger contribution to the 
COBE quadrupole detection, Q rms = 13 ± 6 /uK. In such a case we obtain {a\)r/ {a\)s ~ 



13 ^~ n) and, with the CDM transfer function, b(n) « 0.94 1Q i.285(i-n) ^ ± . 46 ) j 



which, because of the higher cosmic variance, is affected by quite a large error bar. 

Let us consider some examples. From Eq.(6), by taking n = 0.8, we get b xs 2.1; 
even a value of the spectral index quite close to the Harrison-Zel'dovich one, such as 
n = 0.9, involves a remarkable correction by gravitational waves (~ 23%), leading to the 
final estimate b xs 1.3. 

Let us also notice that, given the inflationary model, the COBE data also provide a 
determination of the value of the Hubble constant at the time when the largest observable 
scale left the horizon, i.e. about 60 e-foldings before the end of inflation. For instance, for 
a power-law inflation with p = 11 (n = 0.8), H eo ~ 1.48 x 10~ 4 mp(l ± 0.28), with nip 
the Planck mass. 

An important result of the present analysis is then the possibility to increase the 
estimate of the bias level: this implies less evolution of the considered cosmological models, 
thus lowering the amplitude of pairwise velocities on small scales. The tilted (n < 1) CDM 
models considered here provide a natural solution to the lack of power on large scales and 
excess power on small scales of a CDM model with n = 1 and 6^1. These non-standard 
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CDM models have been analyzed by many authors. Vittorio, Matarrese & Lucchin (1988) 
showed that they imply better agreement with large-scale drifts and with the cluster- 
cluster correlation function. Tormen, Lucchin & Matarrese (1992; see also Tormen et al. 
1992) explored in wider detail the advantages of these models in reproducing the large- 
scale peculiar velocity field traced by optically selected galaxy samples. A good fit of the 
angular correlation function of galaxies in the APM catalog is obtained by Liddle, Lyth & 
Sutherland (1992), with a n ^ 0.5 CDM model. More recently, Adams et al. (1992) have 
considered various cosmological constraints, while Cen et al. (1992) have run numerical 
simulations of n ~ 0.7 CDM models; however, they fix the normalization by fitting the 
COBE data without the tensor-wave contribution. This fact implies an overestimate of 
the small-scale power leading to a residual excess of pairwise velocity dispersion. The 
increase of the COBE determined biasing factor, resulting from our analysis, gives then 
an even stronger support to tilted CDM models. 

While completing this Letter two preprints have circulated that report on independent 
work on similar problems in the frame of various inflationary models (Salopek 1992; Davis 
et al. 1992). By a best-fit of the correlation function, Salopek derives the lower bound 
p ^ 11 (n ^ 0.8) in order to get an acceptable biasing factor. The Davis et al. analysis 
is mostly based on fitting the quadrupole amplitude. Our conclusions are fully consistent 
with their results. 
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Note added The scalar contribution to the multipoles is here calculated through the 
leading term in the Sachs- Wolfe formula, as also done by Smoot et al. in fitting the COBE 
correlation function. However, the next to leading terms (e.g. Abbott & Wise 1984b; 
Fabbri et al. 1987) are no more negligible for £ > 15 and contribute to the smaller angles 
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probed by COBE. Although this effect is not important for the present analysis, it should 
be properly taken into account in fitting the COBE autocorrelation function, especially 
when the experimental sensitivity will be increased. We acknowledge Paul Steinhardt for 
driving our attention to this point. 
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Figure captions 

Figure 1 The ratio of the tensor to the scalar contribution to the CMB multipoles, up to 
order i = 30, as a function of the primordial spectral index n. 

Figure 2 The gravitational-wave correction to the linear bias, b(n)/bo(n) — y/l + D(n), 
as a function of the primordial spectral index n. 
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